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The particle transport through a chain of quantum dots coupled to two bosonic reservoirs is 
studied. For the case of reservoirs of non-interacting bosonic particles, we derive an exact set of 
stochastic differential equations, whose memory kernels and driving noise are characterised entirely 
by the properties of the reservoirs. Going to the Markovian limit an analytically solvable case is 
presented. The effect of interparticle interactions on the transient behaviour of the system, when 
both reservoirs are instantaneously coupled to an empty chain of quantum dots, is approximated by 
a semiclassical method, known as the Truncated Wigner approximation. The steady-state particle 
flow through the chain and the mean particle occupations are explained via the spectral properties 
of the interacting system. 



1 INTRODUCTION 

The advance in technology in the last decades has al- 
lowed the creation of increasingly smaller devices reaching 
the point where the realisation of logic structures on the 
atomic level is possible pj-ljl . Because of their low dimen- 
sionality and temperature the dynamics of the system can 
be dominated by quantum effects, opening a large play- 
ground for experimental testing of many-body correlation 
effects on particle (charge or mass) transport. 

These ideas boosted also the investigation of transport 
of ultracold atoms in systems with reduced dimensional- 
ity. Transport of fcrmionic and bosonic ultracold atoms 
in quantum wires and in one-dimensional optical lattices 
is studied theoretically in [5|-l8[. In [9| a possible re- 
alisation of an atom analogue of an electron quantum 
point contact by the use of a microfabricated magnetic 
waveguide is presented. In an experiment a macroscopic 
atomic cloud was divided into two reservoirs separated by 
a narrow channel by the use of a laser beam [lO| , thus cre- 
ating a cold-atom analog of a mesoscopic conductor. De- 
creasing the dimensionality of the tunneling region from 
one to zero, a new field is investigated - the atomtron- 
ics. The creation of bosonic analogues to the mesoscopic 
systems used in electronic devices like a diode or field- 
effect transistor is suggested in [ll| and also theoretically 



investigated in [12 



In this work we focus on bosonic transport through 
a chain of quantum dots coupled to two bosonic reser- 
voirs that keep the system far from equilibrium. Given 
the by now very well understood behaviour of electronic 
(fcrmionic) systems, the first obvious question is about 
the differences between the bosonic and fcrmionic trans- 
port. It is known that the fcrmionic Anderson impurity 
model - a quantum dot with few energy levels, coupled 
to two electrodes (electron baths) is the simplest possi- 
ble model for a field effect transistor (FET). Since the 
ultracold gas based systems offer a much better degree 
of 'designability' and coherence control, it is also natural 
to investigate the possibility of a bosonic FET. Having 



these goals in mind we offer a formal framework for in- 
vestigation of such systems on the one hand, and on the 
other hand propose a number of efficient and physically 
meaningful approximation techniques, which are able to 
treat even interacting systems. 

In Section 2 we derive a set of stochastic differential 
equations for the time evolution of the reduced system 
by writing down the Keldysh partition function of the 
system and integrating out the reservoir degrees of free- 
dom. In order to derive the set of equations one performs 
essentially the same steps as in [l3| , where a closed system 
is considered, the difference being only in the addition of 
two bosonic reservoirs. In Section 3 we restrict the system 
to the special case of two bosonic Markovian reservoirs, 
which is analytically solvable in the noninteracting case. 
In Section 3.1 we focus on the steady state properties of 
the system. New effects, appearing after an addition of 
an interparticle interaction term to the system Hamilto- 
nian, are explained by the use of the spectral properties 
of the chain of quantum dots. A possible solution in the 
strongly interacting limit is also suggested. In Section 
3.2 the transient behaviour of an initially empty chain of 
quantum dots, which is instantaneously coupled to two 
Markovian reservoirs, is calculated. We find a simple scal- 
ing law between the time needed to reach a steady state 
and the strength of the inteparticle interaction. Section 4 
concludes the paper, offers a possible experimental reali- 
sation of our setup, and outlines the avenues for further 
research. 



2 GENERAL DERIVATION OF A STOCHASTIC 
DIFFERENTIAL EQUATION 

2.1 Single quantum dot coupled to a bosonic 
reservoir 



To start with we consider a system consisting of a single 
quantum dot at energy A coupled to a bosonic reservoir 
with spectral density V{u!) and occupation of the modes 
n{oj). At the initial time U the density matrix of the 



system is assumed to be a direct product of the density 
matrices of the reservoir p and the quantum dot &. The 
reservoir is modeled as a set of noninteracting harmonic 
oscillator levels. Their eigenfrequencies £fc should form a 
continuum, which ensures that the time evolution is irre- 
versible and a steady state is reached. The Hamiltonian 
of the system is given by 

i^ = Aata - E 7fc (S^4 + ^I") + E £fc^I4, (1) 

k k 

where oJtL], create a particle in the quantum dot or in 
the reservoir mode k. One can write down the Keldysh 
partition function [14| . which in the continuous time no- 
tation is given by 

z = ^Y{D[i^*M{LkAt'^)\pk\LkAu)) 



x] D[a.\sL]{a-{U)\a\a+{U)) 

Xe-^k,-iti)-Lk,-.{ti)g-a'_{ti)-a_{ti)^iS^ 

a(i) = {a_{t),a+{t)V. 



(2) 



The —/+ subscript denotes the position of the field on the 
forward/backward branch of the Keldysh contour and the 
ket-vectors \a), \Lk) are eigenvectors of the annihilation 
operators o and Lk- The initial time on both branches 
of the Keldysh contour is denoted by ti and its turning 
point by </. The corresponding action is given by 

S = /,*^dr{at(r)g-i(r)a(r)+ELi.(r)g-i(r)L,(r) 
+ E 7fc (Ll(rKa(r) + at(rKL,(r)) }. 

k 

(3) 
where g ^(r) = {id^ - A)a„ g^ ^(t) = {id^ - e^)a^ and 
a J, is the Pauli z-matrix. 

If one uses the discrete time notation, one can include 
{Lk,-{ti)\pk\LkAti)) e-^^.-(*')-^'=--(*') into the time dis- 



crete form of the matrix g^, "^ (t) 14| and integrate out the 
reservoir degrees of freedom, thus giving the final result 



/^[^ 



,-al(t,)a„(t,) 



\a^{ti)\&\a+{U))e 



iS' 



S' = ///dTldT2at(Tl)G-l(Tl,T2)a(T2). 

G'-^(ri,T2) == S{ti - T2)g-^iTi) - E7fcO-zfffc(n - T2)cr^ 

k 

(4) 
The expectation value of a normally ordered observable 
O = 0{a\a) at the turning point tf of the contour is 
given by 

(ditf)) = jD[a*,a]{{a^{U)\a\a+{U))e--^''y^-^''^ 

xO{aXitf),a_itf)y'''}- 

(5) 
In the same way as in [l3[, where the case of a closed 
system is considered, one can apply the Wigner trans- 
formation {azf{T) = ipir) ± ^ti{t)) and integrate out the 



(6) 



r]{ti),T]{tf) fields to reduce Eq.(l5]) to 

{d{tf)) = J D[rA,ri*,v]{<yw{r{U)A{U)) 

xOwiritf)Aitf))e'^"}. 

The Wigner transform of the density matrix 
<^w{'^*{ti)A{ti)) and the Weyl symbol of the ob- 
servable Ow{ip*{tf),^{tf)) are both obtained after 
integrating out the rf{ti),rj{ti) and ?7*(iy), ?7(iy)-ficlds 
respectively 

(7) 
Calculating Ow is equivalent to writing down the normal 
ordered operator in a symmetrised form and then replac- 
ing a^, a with ip*,ip, respectively. The new action has the 
form; 

S" = I J'f dridT2{77*(Tl)2i(rn + r/2)(Tl - T2)7,(t2) 
+riTimTl-T2){idr, -A) - 2zr(Tl-T2)e(T2-Tl)]77(T2) 

+V*iTiMTi-T2)iidr, -A) + 2ir(ri-r2)e(ri-r2)]VXr2)}. 

(8) 
where r(t) ^ n J^^§^V{uj)j^{uj)e~''^* and {Tn){t) = 

tt/^ |^P(w)7^(w)n(a;)e~*"*. In the noninteracting 
case the action contains only terms that are linear or 
quadratic in the rj, rj* fields. Both types of terms can be 
integrated out to give the following result: 

xjD[rA]{^w{nU)AiU)) 

xoAr{tf)Aitfmhms{hir))}, 

E,,. = 2{rn + T/2){ti-tk). 

. . (^) 
In order to derive the last expression we have divided 

the time interval into N equal parts At = ^ jy ' {ti = 
U + l-At) and defined E G c^+ix^^+i (S;^, = Y.{ti~tk)). 
C,r e C^+i (Ci = ^{ti)). The time evolution of V is 
determined entirely from the argument of the (5-function. 
If one sets /i (ip) to be equal to zero one obtains the fol- 
lowing stochastic differential equation: 

OMt) = -iA^it) - Jl dT2r{t-T)^{T) + c{t), (10) 

where C(i) is a Gaussian stochastic process with zero 
mean and autocorrelation function given by (C(t)C^ [t')) = 
E(i — t'). The equation of motion for ip*{t) is obtained by 
setting /2 (V'* ) equal to zero and it is equal to the complex 
conjugate of Eq. ([TU| . In order to calculate {0{tf)) one 
has to sample a finite number of points {ipjiti)}j=i...NT 
from (TY^!{il;*{ti),ip{ti)), let them evolve according to the 
stochastic differential equation (|10p and then calculate 
the following expectation value: 



Nt 



{o{tf)) « -^T.oMrAf)AAf))- 



(11) 



For large enough tf, a steady state should be reached. 
One should note, that the strength of the memory kernel 
in the second term of Eq. ([TU]) and the autocorrelation 
function of the noise depend entirely on the properties 
of the reservoir. Having a reservoir with constant den- 
sity of states over the entire frequency spectrum, energy 
independent couplings jk and a constant occupation of 
the modes (i.e. T>{lu) = V = const, 7fc = 7 = const, 
ti{lu) = n = const) the Markovian limit is reached, where 
the memory kernel vanishes and the stochastic process 
C^it) becomes a Gaussian white noise: 



dt^{t) = -i^iit)-T^:{t)+C{t) 
{mC{t')) = 2T{n+l)5{t-t') 



(12) 



The addition of an on-site repulsion term ^^a'^a^aa to 
the system Hamiltonian reflects in the action S" by the 
addition of 

(13) 
The r-dependence is dropped for simplicity here. The 
terms in the second bracket are neglected to allow for a 
mapping onto a set of stochastic differential equations. 
This is the essence of the so called Truncated Wigner 
Approximation (TWA) [TJ - flTt . 



2.2 Chain of Af quantum dots coupled to t"wo bosonic 
reservoirs 

The generalisation of the simple example from the 
previous subsection to the case of an arbitrary number 
of wells (quantum dots) A/" between two reservoirs is 
straightforward. The Hamiltonian of this system is given 
by 






J 3 J 



E< 



'k^k^k 



Y^^k^k^k 

k 



Y.lL,ki"'l^k 
k 



Lia,) 



E7fl, 



K^Af) 



E J{ 



;t 



i=i 



"■j+i^'j 



a,,- + aljU 



■j'^j+ij 



kK^Af^k 



j=i 



The ladder operators Lv , R 



(t) m M) 



(14) 
are responsible for 
the annihilation (creation) of an excitation at the left, 
right reservoir and at the chain of quantum dots. We 
always set C/i = = Ufj and Ai = = A^f. 
The corresponding set of stochastic differential equations 
that one has to solve is given by 






~JldT2rLit-T)i,,{T) + iJMt) + CLit) 
~iA^.*^.(t) + zj(^^._i(f) + V^^+i(t)) 



-lU^^p^mnt) 



(1< J < A/") 



dti^j^it) = - /;.dT2rfl(t-T-)VA^(r) + z J^^_i(t) + Cnit) 

(15) 



where TL^fiXL,R are defined in the same way as in Eq. 
pO)) and the subscript L,R refers to the left, right reser- 
voir. We assume that initially the lattice chain is empty 
((ajdj) = 0) and at ij = it is instantaneously coupled 
to the environment, i.e. we take 7j,(i) = lkO{t). The 
Wigner function of the initial state is then 



aw(V'*» = n(fe-2'>^*'^o. 



(16) 



3 RESULTS FOR A CHAIN OF M QUANTUM 

DOTS COUPLED TO TWO MARKOVIAN 

RESERVOIRS 

3.1 Steady state properties of the system 

We first consider the case A/" = 3. Using the nonequilib- 
rium Green's function approach we get exact results for 
the noninteracting case and Markovian reservoirs. The 
mean occupation number Uj [j = 1,2,3) of the dots and 
the steady state current / are given by the following exact 
solutions, for F = 717^2?, A2 = 0, and x = J/F: 



/ 



riL 



TT-L - TiR X" 



2 1 + x2 ' 



TIL + riR 



n2 



nR 



= J- 



1+X2 



nR - ul X 
2 1 + x^ 
{riL - nn) , 



(17) 
(18) 

(19) 
(20) 



where nj^m arc the occupation numbers of the modes of 
the left/right reservoir. We should note that the steady 
state current remains the same independent of the length 
of the lattice chain as long as C/j = = Aj Vj. 

In the case of nonzero interparticle interaction in the 
Markovian limit we approximate the interaction contri- 
bution to the self-energy only by the tadpole diagram 
(one loop diagram with two external legs, also referred 



to as Hartree contribution) [18|. We shall see later that 



already this approximation yields a number of interesting 
details, which are consistent with the TWA predictions. 
In the current case the tadpole diagram renormalises the 
energy level of the middle quantum dot from A2 = to 
A2 = U2{l + nL+ ur) = U2{1 + 2^2). At this point it is 
important to realise that Eq. ([T8| is also valid for A2 ^ 0, 
which means that 77,2 is unchanged in this approximation. 
The same behaviour of 712 is obtained by the TWA. 

From Fig. [Done sees that the current has qualitatively 
the same behaviour in the TWA and in the tadpole 
approximation for not too large 71^. The slope of the 
curves and the position of the peaks in the second ap- 
proximation can be explained with the spectral functions 
of the three quantum dots ^ (a;), (j = 1,2,3) [19| that 
can be obtained from the action S' of the noninteracting 
system after the substitution A2 -^ A2 + 1/2(1 + 2n2)- In 
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FIG. 1: Steady state current for a chain of three quantum dots 
coupled to two Markovian reservoirs for nonzero interparticle 
interaction U2/J = 10^'^. Truncated Wigner approximation 
(triangles - T / J = 5, squares - V / J = 50) and tadpole approx- 
imation (solid line - T/J = 5, dashed line - T/J = 50). Ad- 
ditional parameters: A2/J = 0, ur — 100. The peaks in the 
tadpole approximation are at nz, = 5 x 10"^ and ul ~ 5 x 10* 
for r/J — 5 and F/J — 50, respectively. The corresponding 
new values of A2/J — > A2/J + (72(1 + nL + ue)/ J ~ V^nLJ J 
are 5 and 50. 
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FIG. 2: Spectral functions of the first (left panel) and the 
second (right panel) quantum dot. For F/J = 5 (50) the 
black, dark grey and grey lines in the upper (lower) two figures 
denote the spectral functions for A2/J = 0, 5 and 10 (0, 50 
and 100) respectively. The peak of Ai{y^) is at w = A2. The 
dashed vertical line denotes the value of the critical A2 ~ 
JJin-i where the peak in the current in Fig. [T]in the tadpole 
approximation is reached. 



this approximation the spectral functions of the first and 
third quantum dot are exactly the same since the system 
is symmetric under the exchange of (l,i) <-> (3,-R) 
indices (except ni^_]i) and the retarded Green's functions 
of the system do not depend on rii^^B, in the nonintercat- 
ing case. In the tadpole approximation this symmetry is 
not broken since we have only to renormalise A2. For 



increasing F ^1(0;), ^3(0;) become wider and they do not 
change when varying the energy level A2 of the middle 
quantum dot, except for the appearance of a small dip 
and peak at w = A2. In the following discussion the 
latter effect is not important. On the other hand ^2('^) 
has only a narrow peak at oj = A2. 

Now, we look at the overlap of the spectral functions 
of the left and the middle quantum dots (^i(w),^2(w)) 
(the results for the overlap between A2{^) and .43(0;) 
are exactly the same). For A2 = and increasing F the 
overlap is in the same energy range since the width of 
AiiiS) is almost unchanged in comparison to the width 
of ^1(0;). Only particles in the left dot with energies 
also accessible in A^iyi) can tunnel to the middle dot. 
But this number is smaller since, for larger F, ,41(0;) 
spreads over a wider range of energies and the particles 
at the left dot are distributed over this range. It follows 
that the current should also decrease. This explains the 
difference in the slope of the curves plotted in Fig. [T] for 
small values of n^, — n^. The same behaviour can be 
seen also in Eq. (|20p in the relevant parameter regime 
a; = J/F < 1. 

With this simple picture one can also explain the posi- 
tion of the peaks at the curves plotted in Fig. [1] In the 
tadpole approximation an increase of the interparticle 
interaction strength leads to a change of the energy level 
of the dot (A2 ^ A2 + 1/2(1 + T^L + ur) « U2nL)- We 
have to take into account the competition between two 
effects. On one hand, an increase of n^ leads to a shift 
of the peak of ^2(1^) to higher values, thus decreasing 
the overlap between A\(lij) and An(Lij)^ meaning that 
the relative number of the particles that can tunnel to 
the middle quantum dot decrease. On the other hand, 
looking at Eq. (|17p . the total particle number in the 
first quantum dot increases almost linearly with n^ . The 
position of the peak should be at the point, where the 
first effect begins to dominate over the second one. From 
Fig. [21 we see that for two different F this is the value, 
where ^i(A2 = Vini}i is equal to half of its maximum. 

Within our approximations and keeping the number of 
quantum dots A/" = 3, there is no difference in the results 
for the mean particle occupation Uj (1 < j < Af) of 
the quantum dots in the interacting and noninteracting 
regime. For J\f > 4 such a difference can be seen as shown 
in Fig. [3] for Af = 4 and U2 = U3 = U after applying the 
TWA. The tadpole approximation cannot describe such 
a difference in the particle occupation of the middle two 
dots since it gives the same correction to their energy 
levels A2 and A3. 

We attempted a self-consistent calculation, which leads 
to the following equations for occupations of the middle 
two quantum dots (n2,n3) = (/2(A2, A3), /3(A2, A3)) 
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FIG. 3: Mean particle occupation of the quantum wells for a 
chain of four quantum dots coupled to two Markovian reser- 
voirs, TWA. The black, dashed, dotted and dash-dotted lines 
denote the mean particle occupation in the first, second, third 
and fourth quantum dot. We use the parameters: ul ~ 4000, 
UR = 100, r/J = 5, Ai = Az = A3 = 0, U2 = U3 = U. 



m-. 



/2(f/(n2-l/2),[/(n3-l/2)) = 112 

h{U{n2-l/2),U{n3-l/2)) = n^. 



(21) 



These equations can be solved numerically for a wide set 
of parameters. In the limit of strong interparticle in- 
teractions, we find a better agreement of the emerging 
solutions with the predictions from the TWA for growing 
U. 

To explain the results in the strongly interacting limit 
one has to take into account that each of the Markovian 
reservoirs forces the occupation in the wells to be equal 
to the occupation n^/^j of the reservoir modes. In the 
case TV = 4 and very strong interparticle interactions one 
should expect that the coupling between the middle two 
quantum dots is effectively equal to zero in analogy to 
the self-trapping effect one observes for a Bose-Einstein 
condensate in a double well potential 2l|. One can as- 



sume that the first two quantum dots are coupled only to 
the left reservoir ~ and the last two only to the right one. 
In this case, the occupation of the first two and last two 
dots is equal to nL,nii respectively, which seems to be 
the case after an extrapolation of the results of both ap- 
proximations in the limit of large interparticle interaction 
strengths. 



3.2 Transient behaviour of the system 

In order to find an analytical expression for the be- 
haviour of an empty chain of quantum dots after an in- 
stantaneous coupling with two reservoirs one has to calcu- 
late the retarded, advanced and lesser Green's function 
qR,a,< q£ ^^^ system. The case of a single fermionic 
quantum dot coupled to a reservoir is already considered 



order self-energy (tadpole) approximation. The general- 
isation to a chain of quantum dots and two Markovian 
bosonic reservoirs is straightforward. For U = the re- 
tarded/advanced Green's function is obtained from the 
solution of the set of equations: 

(-za,,-A,)Gf/^(t,i') ^SikS{t-t') + 

(22) 
The retarded/advanced part of the self-energy has the 
form 

^fkit^t') = {-^mt)iSllSkl+6l^f6k^)-J^l^k±l)Sit~t% 

^fkit^t') = {+^mt)iSllSkl+Sl^/6k^)-J^l,k±l)Sit-t'). 

(23) 
After solving Eq. (|22p one can obtain G^ {t, t') by making 
use of the fact that the chain of quantum dots is empty 
at t = 0: 

G<(t,t') ^ /dTidT2G«(t,Ti)I]<(ri,r2)G^(r2,t') 

(24) 
With G^'^'^{t,t') one can obtain all system observ- 
ables. The calculation of the tadpole approximation of 
the Green's functions (denoted by G) of the chain of 
quantum dots is obtained via the following equation; 

Gik^.t') = Gik{t,t') +J:2Uj JjTn,{T)Gij{t,T)G,k{T,t'). 

(25) 
The mean occupation number at the /*'' lattice site is 
then given by 

Mt) = ^D<it,t) 

+ E, 2t/, / drn^ {t)iD< {t, r)Df, (r, t). 

(26) 
The first term is the result from the nonintcracting case 
and the last two are the perturbative corrections from the 
interaction. In the following we consider the case TV = 3 
and observe only the behaviour of n2{t),n2{t). In the 
nonintcracting case we clearly differ between two regimes 
in which the observable has the following form: 



n2{t)^Q.h{nL+nR)fA{t) F < 2^/2 J 



n2{t) = 0.5(ni + ni^)fB{t) F > 2^/^J 

with fA{t),fB{t) given by: 

fA{t) = l + ^(-8j2+F2cos(t/3)-F/3sin(t/3)) 
fsit) = l + ^(8j2-F2cosh(t/3)-F/3sinh(i/3)) 



(27) 



/? = ^|8J2-F2| 



(28) 



In the regime T > 2^/^ J (Fig. |4]) the observable 
converges exponentially to its steady state, as in the 
case for the particle occupation of a single quantum dot 
coupled to a Markovian reservoir. The time scale of this 
process is proportional to r/(4J^). But in the limit of 
very small F, one observes a step-like behaviour of the 
particle occupation, the length of the steps being equal 
to 2tt/(3. One can also see that the fastest convergence 
to a steady state is obtained in the case where T ^ J. 

The next task is to see if the interparticle interactions 
at the middle dot can influence this transient behaviour. 
For the special case of TV = 3 one can bring Eq. (P^ into 
the more compact form 

Mt) = n2{t)+4U2jdTn2{T)^G^2{t,T)tG<^{T,t)). 

(29) 
The correction to the particle occupation in the mid- 
dle quantum dot is zero. It follows that not only the 
steady state but also the transient behaviour of n2(t) 
is unchanged by the presence of interactions at least 
within these approximations. The situation is differ- 
ent if one looks at the numerical solution of Eq. (fT5|) . 
where all classical contributions of the interparticle inter- 
action are taken into account. In both parameter regimes 
(F ^ 2^/^ J) one observes a quadratic dependence of the 
time needed to reach a steady state from the interparticle 
interaction in the middle quantum dot (Fig. |4]). 
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4 CONCLUSIONS 

Wc have studied the transient behaviour and the steady 
state properties of a chain of quantum dots that is instan- 
taneously coupled to two Markovian reservoirs. For the 
case of three dots an exact solution in the nonintcract- 
ing case is shown. We see that the interparticle interac- 
tion does not change the mean particle occupation in the 
middle well in both the TWA and the tadpole approxima- 
tion. But the time the system needs to reach a stationary 
state increases quadratically with the interaction in the 
TWA. We have also found a qualitative explanation for 
the behaviour of the steady state current by the use of 
the spectral properties of the chain of dots. Increasing 
the number of wells from three to four, additional effects 
arise from the interparticle interactions. Here the interac- 
tion effectively reduces the coupling between the middle 
two dots such that n\ = 77.2 = n^ and nj, = 774 = riji in 
the limit of very strong interactions. 

In order to access this interesting physics experimen- 
tally we envisage the following procedure, which has es- 
sentially been partly realized already by the authors of 
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One starts with a rather large trap with a Bosc- 
Einstein condensate in perfect equilibrium in it. Then 
by an instantaneous potential shift one induces a slosh- 
ing of the condensate. After that the system shoud be 
cut into two subsystems, for instance by an impenetrable 



FIG. 4: Time evolution of /A(i),/s(t) defined in Eq. ((28]) 
for [/ = 0, r = ^2'^" J (solid line in the upper panel) and 
r = 5J (solid line in the lower panel). The other lines are the 
results from the TWA obtained after dividing n-iif) by (til -\r 



nR)/2. The values of U2/J are 5 x 10 



,10- 



,10"^ (10^^ 5 X 



lO^'*, 10^^) for the dashed, dotted and dashed-dotted lines in 
the upper (lower) panel. In the inset one can see the time that 
fA(t) or fB{t) needs to reach 0.95. The numerical results are 
fitted with a curve of the form g{U2) = a + bUi- 



barrier. In this way one produces two different bosonic 
reservoirs which contain a large number of particles in 
excited states. Gradually removing the barrier one can 
then couple these "reservoirs" and hence allow for the 
transport. The additional structuring of the contact area 
into several quantum dots can be accomplished in the 
way similar to that described in |ll| for one well, or by 
adding a lattice potential along the channel created in 
[10|. Wc hence expect that such a 'bosonic FET' can 
be manufactured with the state-of-the-art experimental 
methods. 

Needless to say, there is enough room for improvement 
of our approach. While an extension of the TWA ap- 
pears to be highly non-trivial, the inclusion of the higher 
order self-energies is, in principle, rather straightforward. 
Since the latter will definitely generate energy-dependent 
quantities, we expect not only quantitative but also qual- 



itativc differences to our predictions to emerge. However, 
they would only play a significant role for intermediate 
to strong interactions. 
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